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We propose a minimal model for the Josephson current through a quantum dot in a Kondo 
regime. We start with the model that consists of an Anderson impurity connected to two 
superconducting (SC) leads with the gaps A^ = |Aa| e'"", where a = L, R for the lead at left 
and right. We show that, when one of the SC gaps is much larger than the others |Al| 3> | A_r|, 
the starting model can be mapped exactly onto the single-channel model, which consists of the 
right lead of A_r and the Anderson impurity with an extra onsite SC gap of Ad = Fl e' ^ . 
Here 6l and Tl are defined with respect to the starting model, and Tl is the level width due 
to the coupling with the left lead. Based on this simplified model, we study the ground-state 
properties for the asymmetric gap, |Al| ^ I^hI, using the numerical renormalization group 
(NRG) method. The results show that the phase difference of the SC gaps ({> = 9r — 6l, which 
induces the Josephson current, disturbs the screening of the local moment to destabilize the 
singlet ground state typical of the Kondo system. It can also drive the quantum phase transition 
to a magnetic doublet ground state, and at the critical point the Josephson current shows a 
discontinuous change. The asymmetry of the two SC gaps causes a re-entrant magnetic phase, 
in which the in-gap bound state lies close to the Fermi level. 

KEYWORDS: Kondo effect, Josephson effect, quantum phase transition, Anderson model, numerical renor- 
malization group, quantum dot 



> 

oo 

o 
o 



o 



X 



1. Introduction 

The Kondo effect in superconducting (SC) materials 
was originally studied for dilute magnetic alloys. ^"^ In 
normal metals, the ground state of the magnetic impu- 
rity is known to be a spin singlet, as a result of the an- 
tiferromagnetic coupling with the conduction electrons.** 
However, in the presence of the SC long-range order, the 
energy gap of the host material disturbs the conduction 
electron screening of the local moment of the magnetic 
impurities. It was shown in early years that whether the 
local moment is screened or not is determined by the 
ratio of the Kondo temperature Tr- to the SC gap A. 
The ground state becomes a doublet for A S> Tk, but it 
is still a singlet for A <C Tk- In contrast to the normal 
metals, for which the magnetic solution of the mean-field 
theory^ does not describe the correct ground state, in the 
superconducting case the bound state appearing in the 
energy gap causes the occurrence of a quantum phase 
transition to a magnetic doublet ground state. 

These aspects of the Kondo physics in superconduc- 
tors have been re-examined precisely by efficient numer- 
ical methods such as the quantum Monte Carlo (QMC)® 
and Wilson numerical renormalization group (NRG) ap- 
proaches. ^"*° Furthermore, it has also been reconsidered 
for novel systems such as quantum dots.^*'*^ One of the 
new features in the quantum dot coupled to two SC 
leads is that the ground-state properties of the system 
can be controlled by the phase difference of the two SC 
gaps (j). It induces the Josephson current flowing through 
the dot, and at the same time it affects the screening of 
the local moment. When the local moment remains un- 
screened, the spin-flip tunnelling through the magnetic 



impurity causes a current flowing in the opposite direc- 
tion to that in the case of the normal junctions. An ex- 
planation has been given based on the lowest-order per- 
turbation theory with respect to the spin-ffip tunnelling 
Hamiltonian or equivalently treating the magnetic mo- 
ment to be a classical spin.*^"^^ However, in the Kondo 
regime, the quantum-mechanical nature of the moment 
is essential to the screening. Therefore, in order to clar- 
ify how the Josephson current and the dynamics of the 
moment affect each other, the higher order terms in the 
tunnelling matrix element must be taken into account. So 
far, several studies have been carried out using the non- 
crossing approximation (NCA),^^'*^ slave-boson mean- 
field (SBMF) theory,!^'!^ perturbation theory in [7,^° 
QMC,^* and NRG.^^ Rozhkov and Arovas compared the 
ground-state energy for (jj — and (jj — tt to discuss 
the phase transition from the to tt phase. *^ Clerk and 
Ambegaoka pointed out that this transition can be ex- 
plained in terms of the bound state appearing in the SC 
gap,*'' as in the case of the magnetic impurity in the 
bulk superconductors.*"'^ So far, mainly one special case, 
in which the absolute value of the two SC gaps are the 
same jA^j = \Aji\ and the couplings to the two leads 
are equal Fl — F^, has been examined. *^"^"'^^ However, 
the asymmetry in the two gaps \Al\ ^ \Ar\ and that in 
the mixing parameters F^ ^ F/j affect sensitively the SC 
proximity effects on the impurity site. The superconduct- 
ing correlations penetrating from the two leads determine 
the effective field that is coupled directly to the impurity. 
Therefore, the asymmetries in the gaps and in the mixing 
matrix elements, as well as the phase difference </>, will 
also affect the bound state and low-temperature proper- 
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tics of the impurity. 

The purpose of this work is to study the interplay 
of the Josephson current and the dynamics of the lo- 
cal moment. To this end, we start with the Anderson 
impurity that is connected to two SC leads. Although it 
is already a simplified model, it still contains a number 
of parameters. In the present paper we mainly consider 
the asymmetric gap of jA^j ^ \^r\, where one of the 
SC gaps is much larger than the others. We show that 
in the limit of \Ai\ -^ oo the starting model can be 
mapped onto a single-channel model, in which the left 
lead is chopped off leaving an extra onsite SC gap of 
Arf = r^ e '^ at the impurity site. This single-channel 
model possesses both the Josephson and Kondo aspects 
since the phase difference can be defined with respect to 
Arf and A^, and the simplification gives us a great ad- 
vantage in carrying out numerical calculations. Based on 
this single-channel model, we examine the ground-state 
properties for jA^j » \Ar\ using the NRG method. The 
phase diagram of the ground state is calculated for a 
wide parameter range of the Coulomb interaction U and 
the mixing strength F^. The results show that the phase 
difference 4> tends to make the stability of the singlet 
ground state worse, and it can hold a casting vote to de- 
termine whether the ground state is a singlet or doublet 
when the other parameters are competing. Specifically, 
at ~ TT, r^ ~ Ffl, and near half-filling, the bound state 
lies close to the Fermi level, and it causes a re-entrant 
behavior of the magnetic doublet ground state. 

In §2, the starting model and the mapping onto the 
single-channel model are described. In §3, the bound 
state in the noninteracting case is discussed briefly. In §4 
the results of NRG calculations arc presented. Summary 
is given in §5. In Appendix, the perturbation theory for 
a zero-energy bound state and other related matters are 
provided. 

2. Mapping onto a single-channel Model 

In this section we introduce the single-channel model, 
which captures the essence of both the Josephson and 
Kondo physics, starting from the model consisting of a 
single impurity and two SC leads. 

The Hamiltonian of the Anderson impurity connected 
to two superconducting leads, at left (L) and right (R), 
is given by 



n = n,! 



njU 



2_, ( Wa + 'H^x 



(1) 



\=L,R 



where TiJ] and Ti^ represent the impurity part, and A = 
L, R is the label assigned for the SC leads. The explicit 
form of each part is given by 

n'd = U + ^)ind-l) , n^^jin,-lf, (2) 



'^^ — 7 .^Afc^Afcg'^AfcCT 



ka 



2^ ( ^AC^fcl C^_fei + ^ACA-fei<^AfcT 



), (3) 



H^ 



dX 



E^> 



di 



(4) 



The operator dj. creates an electron with spin a at the 
dot, and n^ = ^^dl.d^. Similarly, c^^,^ creates a con- 
duction electron in the lead at A, and e\k and A\ are 
the energy and SC gap for the electrons in the lead, re- 
spectively. The tunnelling Hamiltonian 7Y J_^ connects the 
dot at the center and lead at A. At the interface the lo- 
cal operator of the conduction electrons is defined by 



'^Acr ~ X]fe '^Afccr/ V A^- Correspondingly, the current fiow- 
ing into the dot from the left lead Jl, and the current 
flowing out to the right lead Jn, are given by 



Jl ^ iev^'^i 4 



-^cr'-La ^La"'cr 



Jr = '^eVj^Y^ 



"Ra^'cr "'a^Ra 



(5) 



(6) 



Note that we are using units fi = 1 unless otherwise 
noted. 

To discuss the superconducting properties, we use 2x2 
matrix Green's function of the Nambu formulation. 



Gdd{T) 



{TrdjT)dh {TrdAT)d{) 
{Trd\{T)d\) {Trd\{T)d{) 



(7) 



The Fourier transform of the impurity Green's function, 
Gddi^^n) — /gdre"^""^ Gdd{T), can be expressed as 

{Gd<i(i^«)}" = it^nl - EdTz - vlTzgL(iuJn)Tz 

~ «flT3fl'fl(i^«)T3 - T.{\UJn) , (8) 

where Ed^ €d + U /2, and Ti for i = 1, 2, 3 is the Pauli 
Matrix, and S(iw„) is the self-energy due to TL^. Also, 
the Green's function for the bulk superconductor is given 
by {g\k{^^n)y = it^nl - CAfcTa - Aa. The interface 
Green's function which enters in Eq. (8) is defined by 
g.^{iujn) = X]fefl'Afe(i'^n)/^i s-iid can be written as 

9 a (i^") = -T^PxiO) 



V^4+W^' 



(9) 



where the density of states Pxi^) — J2k ^(^ ~ ^\k)/N has 
been assumed to be a constant. The SC order parameter 
enters in the off-diagonal element of Aa, 





At 



Aa 




= |Aa| (cos^aTi -sin6'AT2' 



(10) 
where A a — |Aa|c'^^. When the phase difference = 
Or. — Ol is finite, the Josephson current flows through the 
impurity. Furthermore, due to the hybridization with the 
electrons in the SC leads, the superconducting correla- 
tion Xd = (did-^), which corresponds to the off-diagonal 
element of Gdd, evolves at the impurity site. 

In the noninteracting case, Ti,^ 
Green's function is given by 
-1 



0, the impurity 



{GS,(i^„)}" = 



iuj„l- EdTs 



iujnl - A J 



iw„l - Ar 



v^^^ 



^Rr 
(11) 
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where T^ ~ Trp^(0) vj^. Correspondingly, the onsite SC 
correlation and Josephson current can be written as 

1 



x2 



-T — 



TrA, 



TbA 



R^R 



{Jr)" 



-1 



(12) 



-T — 

/3trdet{G°,(ic.„)} 



4r^ri|A, 



^l\ sm[Oji 



-Ol) 



Vc.2 + |A^IV^U|Al|^ 



(13) 



where /3 = 1/T, and (Jl)" = {Jr)" . The Coulomb 
interaction U disturbs the SC correlation to penetrate 
into the impurity. This tendency can be seen already in 
an approximation of the Hartree-Fock level, through an 
anomalous contribution^'^ 



n_/U 



U 



7t^t 



twt\ 



d\dl (d^d^) + {d^^d\)d^d 



(14) 



This term corresponds to a Fock term in the Nanibu 
formulation 



U 





X*d 



Xd 




(15) 



and has a different sign from that of the attractive in- 
teraction which drives the leads superconducting states. 
However, to study precisely the effects of the Coulomb 
interaction on the Josephson current and on the SC cor- 
relations at the impurity, the higher-order terms beyond 
the mean- field theory should be taken into account. To 
carry out nonperturbative calculations, the present ver- 
sion of Anderson impurity already contains a number of 
parameters compared to that in the normal leads, i.e., 
we have |A;^| and 9\ in addition to T\, e^ and U. There- 
fore, it seems to be meaningful to examine some special 
cases, at which the model can be simplified without los- 
ing the essence of Josephson and Kondo physics. In the 
following, we consider some such special cases. 

2.1 Model I: \Al\ ^ \Ar\ 

One of the cases that has been studied by a number of 
authors so far is |Ai,| = \Ar\ (= A), where the absolute 
value of the two gaps are equal. Physically, it means that 
the two leads are made of the same SC material. In this 
case the noninteracting Green's function can be written 
in the form 



{Gl5rf(ia;„)} ' =iLUnl-EdT3 + {rL+TR[ 



(16) 



A 



Tl+Tr 



{Tl COS0L + Tr COS0r) n 



(r 



L sm^L 



r 



R SmtiR) T2 



(17) 



and thus we have only one characteristic energy scale 
for the superconductivity, i.e., A. However, in this case 
the two-channel nature of the system is still preserved 
for finite 6. This feature can be seen in the functional 



form of G'^^{iujn) in Eq. (16). It becomes equivalent to 
that of a single-channel model only for e*^'* = e*^^ , when 

~ 2 

the amplitude of A coincides with A, i.e., A = A^ 1. 
Therefore, in the case of jA^j = \Ar\, both of the su- 
perconductors must be taken into account explicitly to 
investigate the Josephson effect. 

2.2 Model II: \Al\ ^ \Ar\ 

The situation is different when one of the gaps is much 
larger than the other's, |Al| ^ \Ar\. Specifically, in the 
limit of \Al\ -^ oo, the noninteracting Green's function 
G^^iiujn) can be written as 

{G°,(ic.„)}-' = 

iw„l - EdT3 - TlTi - v\ X3 QRilUJn) T3 . (18) 

where the phase of A^ is chosen to be 61^ = 0. One no- 
table feature is that the third term in the right-hand side 
of Eq. (18) can be regarded as a static superconducting 
gap induced at the impurity site, and its value is given by 
r^ . This term comes from the mixing self-energy due the 
left lead v\T^gj^{iujn)T^, and in the limit of \Al\ -^ oo 
the off-diagonal SC correlation remains finite to pene- 
trate into the impurity site while the diagonal damping 
part corresponding to the level width vanishes. There- 
fore, in the case of \Al\ ^ \Ar\ the information about 
the left lead can be included through the off-diagonal 
FlTi term of the impurity Green's function, which can 
also be described by an effective single-channel Hamilto- 
nian with an extra SC gap at the impurity site A^; = F^ , 



nj 'T_/0 

/tcff — rid 



o/SC 



= Ad 



njU 

- n,d 

d\d\ 



i/SC 

n-d 



d^d^ 



Hr + H, 



dR 1 



(19) 
(20) 



The phase difference (f> defined with respect to the start- 
ing model is described by the phase difference between 
Ar and Ad. Using the effective Hamiltonian Ticff, one 
can calculate the expectation values such as the Joseph- 
son current (Jr) and the Green's function Gddi^^n) for 
the interacting electrons, which coincide with that of the 
starting model in the limit of \Al\ — s- oo. The map- 
ping introduced here is exact. It can be confirmed, for 
instance, using the path integral formulation: the same 
effective action for the impurity can be derived from both 
Ti. and Tics after carrying out the integration over the 
conduction-electron degrees of freedom [see Appendix 
A] . The single-channel Hamiltonian Tioff can be regarded 
as a minimal model that possesses both the Josephson 
and Kondo aspects of the system. This simplification 
gives us a great advantage in carrying out numerical cal- 
culations. Particularly, the NRG approach works excel- 
lently for single-channel systems, while the numerical ac- 
curacy becomes rather worse for multi-channel systems 
because the number of the low-lying states that should 
be retained in the calculation increases with the num- 
ber of the channels. Therefore, in the present study, we 
mainly consider the asymmetric SC gap of |Ai| 3> \Ar\ 
using the single-channel Hamiltonian Eq. (19). 
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Fig. 1. Phase diagram of the ground state for (p = n, Tj^ = Fr 
(= F), and E^ = 0. Here \Al\ is chosen to be |Ail > |Afl|. 
An infinitesimal U lifts the 4-fold degeneracy caused by the zero 
mode to stabilize the singlet or doublet state depending on the 
value of lAfll/IA^I and l/F. The results are obtained with the 
perturbation theory in Appendix D. 



3. Bound state in the gap 

The bound state appearing in the SC gap plays an 
important role on the ground-state properties of dilute 
magnetic alloys. ^"'^ In the case of the quantum dots, the 
Josephson phase (/) — 9r~9l also changes the energy and 
wavefunction of the bound state, "'^^ and it also affects the 
screening of the local moment. To make these features 
of the (j) dependence of the bound state clear, we discuss 
briefly about the bound state in the noninteracting case 
[see also Appendix B]. 

The bound state can exist in the gap region |e| < 
mindA^I, |Ai^:|), and the eigenvalue is determined by the 
equation det|G2^(e)} — 0. For noninteracting elec- 
trons the determinant can be written explicitly as, 



Aet{GU^)Y' ^2[F2{^) - F,{^)] 



\{^') ^ 


^[rl+T% + E^,-e^] , 


ofe^l - 


Tl^' Tne^ 




VIAlP-e^ y/\An\^-e^ 




TLTnie^ - |Al||A^| cos</)) 



(21) 



(22) 



V|A,.P-,2^A^2 



Here i^i(e^) is a simple decreasing function of e^, while 
the functional form of F2{e^) depends on the parameters 
jAj^l, |A^| and 0. However, at e = 0, it is given simply 
by F2(0) = -TlTr cos(/). Therefore Fi{0) > i^2(0), and 
the equality holds for cfi ^ tt, Tl = Tr and Ed — 0. We 
consider this particular case in the following, and then 
discuss briefly the bound state for the models mentioned 
in §2. 

3.1 Zero mode for (j) — n, Tr = F^ and E^ = 

We consider here the special case where the three con- 
ditions hold simultaneously; i) vr-junction (j) = n, ii) 
equal coupling Tl = F^ (= F), and electron- hole sym- 



metry Ed ~ 0. In this case -P^le^) is a increasing function 
of e^, and at e = it coincides with -F'i(O), so that the 
bound state comes just on the Fermi energy e = 0. 

Thus, the ground state has 4-fold degeneracy with re- 
spect to the occupation of the zero-energy bound state 
by the Bogoliubov quasiparticles [see also Appendix C] . 
The Coulomb interaction H^ lifts the degeneracy. We 
have calculated the energy shift for an infinitesimal posi- 
tive U using the perturbation theory, the outline of which 
is given in Appendix D. The ground state is a singlet 
for F > Fcr, while it is a doublet for F < Fcr. The critical 
value Fcr depends on the ratio x = |Aii;|/|Ai,| as shown 
in Fig. 1. In this figure the phase boundary is plotted 
in the x vs l/F plane choosing |Ai| to be |Al| > |Aij|, 
where the energy is scaled by |Ai,;|. The critical value 
Fcr increases with x and diverges exponentially at x = 1, 
and thus for |Ai| = |Afl| the ground state is a doublet 
independent of the value of F. This is because in this 
particular case the SC correlation penetrates from the 
left lead and that from the right cancel each other out to 
make the net SC correlation at the impurity site Xd to be 
zero. This example shows that the asymmetry in the two 
gaps [Ail j^ \Ar\ affects the screening of the moment 
crucially when the other parameters are set to be highly 
symmetric. Note that in the perturbation theory an ex- 
change integral (Fock term) of the form Eq. (12) shifts 
the bound state from the Fermi level and favors the sin- 
glet state, whereas a direct type contribution (Hartree 
term) favors the doublet state. 

3.2 Bound state for Models I & II 

When the absolute value of the two gaps are the same, 
[Ail = [AflJ = A, Eq. (23) simplifies as 



(Fl + Tr) e^ 
VA2 - £2 



-FrF 



Lifl 



2FiFj^A2 sin^l 
A2-e2 ■ 



V l^^ — £" AA- — fc- 

(24) 

Since this is an increasing function of e^ which diverges 
positively at e^ -^ A^, the bound state exists for any 
finite A. 

For the other model with the asymmetric gap |Ai| ^ 
|A/j|, the function i^2(e^) is given by 
/I 



F2ie') = F, 



- Arf|Afl| cos( 



(25) 



(23) where A^ = F^. In this case the behavior of i^2(e^) at 
e^ -^ I A-rI^ — 0^ depends crucially on the sign of the nu- 
merator, and thus whether or not the bound state exists 
depends on the values of the parameters. 



4. NRG approach 

In the rest of this paper we discuss the ground state 
properties of the Anderson impurity in the limit of 
[Ail ^ |Aij|, where one of the SC gaps is much larger 
than that of the others, based on the single-channel 
model given by Eq. (19). As mentioned in §2.2, this 
model contains essential aspects of both the Josephson 
and Kondo effects. Also, the reduction in the number of 
the channels gives us a practical advantage in the NRG 
calculations. 
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4.1 Method 

Through a standard procedure of the logarithmic dis- 
cretization,^"^ the conduction band can be modeUed by 
a hnear chain with the complex pair potential A^ ~ 
I AjjI e"^, and a sequence of the NRG Hamiltonian H^ is 
given by 

Hn =A(~-i)/2 [7^0 + H^+ nT + H,R + Hr] , (26) 

1 + 1/A 



Hn=D 



N-l 



7 , 7 , S" ^ ( Jn+l a ilia + J na Jn+1 a j 

1=0 a 

N 
E {^Rfyli+^*Rfnifn^) > (27) 



n=0 



HaR=DY,^R{fld.+dlh. 



(28) 



Here D is the half-width of the conduction band, and 
fna is an operator for electrons in the right lead. The 
hopping matrix elements ^„ and Vj^ are defined by 

1 - l/A"+i 



£,n = 



v/1 - 1/A2"+Vl - 1/A^"+^ ' 



2ri^AA 



A, 



1 1 + 1/A 



(29) 



log A, (30) 



ttD ' "'' 2 1 - 1/A 

where Aa — > 1, in the continuum limit A — > 1.23,24 ^]^q 
low-lying energy states of Heff can be deduced from that 
of A^'^^^^/'^iJjv for large N . In the presence of the com- 
plex SC gap Hm does not conserve the charge, and the 
total spin 5* is the only one quantum number that can 
be used for the block diagonalization. In the two special 
cases, = and tt, the Hamiltonian Hjq has an extra 
U(l) symmetry corresponding to the x-th component of 
the axial charge I, = E1-i(-1)" {f^fU + fnif^)/'^ 
with f_i^ = dj.,^'^ and the SC lead can be transformed 
into a normal lead with a staggered potential J 

The ground-state average of the Josephson current can 
be obtained using the discretized version of Eq. (6), 

{J)n = 



,DY, U^N\lld,\^N) - ($iv|4/oJ* 



N, 



(31) 



where |$Af) is the ground state of iJjv- The expectation 
value can be calculated successively for large N using the 
recursive relation among the matrix elements {J)n and 
(J)iv+i- We have carried out the calculations retaining 
the lowest 500 states and taking A and \^r\ to be A = 2.0 
and \1^r\/D = 1.0 x 10~^. In the following, we use |Aii;| 
as the unit of the energy in most of the figures, since that 
is a typical energy scale for the superconductivity. 

4-2 Results at half-filling 

In this subsection we show the results obtained in the 
electron-hole symmetric case, e^ — —U/2. We first of all 
consider the ground state for A^ — 0. Since the value 



u 




5 10 15 20 25 30 

U/A 



Fig. 2. Phase diagram of the ground state for the symmetric An- 
derson model e^ = —U/2 connected to the single SC lead of A, 
where the SC gap in the impurity site is taken to be A^j = 0. 
The NRG calculations have been carried out for A = 2.0 and 
A/D = 1.0 X lO-'^. 



of Ac( is defined to be equal to T^ in the original two- 
channel model Eq. (1), the condition Arf = means that 
the left superconductor is disconnected completely from 
the impurity. In Fig. 2, the phase diagram of the ground 
state is shown as a function of C//|A/j| and rij/jAflj. The 
ground state is a spin-singlet state for large T r or small 
C/, and it is a doublet for small Vr or large U . We note 
that this particular case, A^ = 0, has been examined by 
Yoshioka and Ohashi,^ and the phase boundary is deter- 
mined basically by a single parameter Tk/\^r\ in the 
Kondo regime. Nevertheless, for small U or the electron- 
hole asymmetric case, in which the contributions of the 
charge excitations still remain, the low-energy proper- 
ties are characterized not only by the single parame- 
ter T/f, but the Wilson ratio R and the renormalized 
impurity level e^ are necessary for a complete descrip- 
tion.*'^^ Therefore, in order to see an overall picture of 
the ground-state properties, the phase diagrams which 
are plotted as a function of the bare parameters, as Fig. 
2, are convenient. 

When the SC gap of the impurity /S.d (= F^) be- 
comes finite owing to the coupling with the left lead of 
|Al| — !■ 00, the phase difference cj) can be introduced 
at the interface of A^ and A^. We now examine how 
(j) affects the ground state in the equal-coupling case 
Fl — Vr- In Fig. 3 the phase diagram of the ground 
state is plotted as a function of U and F^ for several 
value of the phase difference = 0, 7r/2, 37r/4, tt. Here 
Arf is varied simultaneously with T r keeping the relation 
Ad = Vr. We see that tends to disturb the screening 
of the local moment and enlarges the doublet phase in 
the figure. The feature of the phase diagram at 7^ tt is 
similar to that for A^ = in Fig. 2. However, the fea- 
ture is quite different for = tt. This is caused by the 
bound state lying just on the Fermi level at fj = 0. An 
infinitesimal U lifts the degeneracy to stabilize the sin- 
glet state at Tr > Fa, and doublet state at Tr < Fcr- In 
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Fig. 3. Phase diagram of the ground state for the equal coupUng 
r^ = Ffl at half-fiUing e^j = — {//2 is plotted for several values 
of the Josephson phase ip. Here A^ (= F^) is taken to be equal 
to Ffl. For = TT and f7 = 0, the phase boundary is given by 
Fcr = 1.77064|Afl|. This special behavior at </< = tt is caused 
by the zero mode, and an infinitesimal U stabilizes the singlet 
(doublet) ground state for Fn > Fcr (Fij < Fcr)- 



the present model, A^ — > cx), the critical value is given 
by Tcr = 1.77064 1 Afl I as described in Appendix D. For 
this reason, in Fig. 3 the phase boundary for cf) — n starts 
from the finite value, Fcr, in the y axis. 

To see how the singlet or doublet ground states evolves 
in the successive NRG steps, in Fig. 4 the low-lying en- 
ergy levels of A~(^~^'/^iJ7v is plotted as a function of 
N for (a) U = \Ar\, and (b) U == 9.0|Aii|. Here the fac- 
tor A"*^^"-"^'/^ restores the original energy scale of Tieff, 
and the ground-state energy is subtracted from the eigen- 
value for the excited states. The parameters are taken to 
he (j) = n and Tr = 1.5\Afi\. Specifically, the SC gap 
in the impurity site is chosen to be A^ — \Aj^\, so that 
the couplings are asymmetric F^ ^ Tji in this parameter 
set. In the figure the eigenstates are labeled by the two 
quantum numbers {Ix, 25), where Ix is the U(l)-axial 
charge mentioned in the above and S is the total spin. 
The solid and open symbols correspond to the levels for 
even and odd N, respectively. For large N, the excitation 
energies converge to the fixed point values, and have no 
even-odd oscillatory dependence on N, which is a typi- 
cal behavior seen in the presence of the energy gap.^ The 
ground state is a singlet for A^ > 32 in (a), while it is a 
doublet in the case of (b) . In the figures all the low- lying 
levels are shown for E < |Afl|, but not all are shown 
for E > |Afl|. The levels below the gap |Afl| can be 
classified according to the occupation of the bound state 
with the Bogoliubov quasiparticles, which can be inferred 
from the value of Ix because it is transformed into the 
usual charge by the Unitary transformation described in 
Ref. 7. In the case of (a), where U is relatively small, the 
first and second excited states lie in the gap region. The 
occupation of the bound state in the ground, first, and 
second excited states are supposed to be empty, single, 
and double, respectively. In contrast, in (b) the occupa- 
tion of the bound state in the ground and first excited 
states are supposed to be single and empty, respectively. 
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Fig. 4. Low-lying energy levels of A~^^~^'''^Hpf as a function 
of N, where the energy is scaled by jA^I and is measured from 
the ground-state energy at each N. The solid and open symbols 
are the levels for even and odd N, respectively. Only the lowest 
eigenvalue in each subspace labeled by the quantum numbers 
{Ix, 2S) are shown, where Ix is the U(l)-axial charge and S is 
the total spin. The parameters are taken to be t/i = 7r, A^ = 
l^fll. ^R = l-5|Afl|, ea = -U/2, and A = 2.0. At large N, the 
levels approach to the fixed point values showing no even-odd 
dependence on N. For A'^ > 32, the ground state is a singlet in 
(a) U = |Ah|, while it is a doublet in (b) U = 9|Afl|. 



The doubly occupied state could not stay below the gap 
in (b) because of large U. 

The eigenvalues depend on the phase difference (j), and 
particularly the fixed-point energy at large N determines 
the low-energy properties. In Fig. 5, the fixed-point value 
of the lowest excitation energy Es=o — Es=i/2 is plotted 
as a function of U for several values of (/>, where Es=o and 
Es=i/2 are the lowest eigenvalues of A"'^^"-'^)/^ Hn in the 
singlet and doublet subspace, respectively. The parame- 
ters are chosen to be F_r — 3.77697|A_r|, and Ad = A^. 
For U < Uc with Uc — 25|A_r|, the ground state is a 
singlet, and the excitation energy depends visibly on (p. 
The energy separation between the two states decreases 
with increasing 0, and it means that (j) tends to desta- 
bilize the singlet ground state. On the other hand the 
ground state is a doublet for U > Uc, and the (j) depen- 
dence of the excitation energy is very weak. This seems 
to be caused by the fact that the local moment, which 
remains in the doublet ground state, makes the SC cor- 
relation length small, and the phase coherence between 
the dot and SC lead becomes weak. In the singlet state, 
however, the phase coherence is preserved, and thus the 
excitation energy depends sensitively on (j). We can also 
see in Fig. 5 that the critical value Uc decreases with in- 
creasing (f) to reduce the area of the singlet ground state 
in the phase diagram. 

For asymmetric couplings Tr ^ F^ , the phase diagram 
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Fig. 5. The fixed point value of the lowest excitation energy of 
\-(N-l)/2fj^ at large Af as a function of U. Here -Bs=0 a^nd 
Eg— 1/2 are the eigenvalue of the lowest singlet and double states, 
respectively. Therefore, the ground state is a singlet (doublet) for 
negative (positive) Es=0 ~ ^S=i/2- The parameters are taken to 
be Aa = \An\, Tr = 3.77697|Afl|, ea = -U/2, and A = 2.0. 



of the ground state has some different features from that 
of the equal coupling case seen in Fig. 3. Two typical 
examples for the asymmetric coupling are shown in Fig. 
6 for several values of </>. Here the value of A^; is fixed 
to be (a) |Afl| and (b) 2|Aj^| independent of the cou- 
pling Ffl. In each of the figures (a) and (b), all the lines 
start from the same point on the horizontal line, the po- 
sition of which is given hy U — 2Ac; from the atomic- 
limit solution at F^ = 0. A re-entrant behavior seen for 
~ TT is caused by the bound state lying near the Fermi 
level. The phase boundary for cj) ~ t: crosses the vertical 
axis at F^ = A^, at which the couplings becomes equal 
Fi — Tji and the bound state comes just on the Fermi 
level. An infinitesimal U lifts the 4-fold degeneracy of 
the ground state as shown in Fig. 1, and specifically in 
the present case, A^ — > oo, the critical value is given by 
Fcr = 1.77064|Afl| (see Appendix D). Thus the exam- 
ples (a) and (b) represent typical cases of A^ < Fcr and 
Ad > Fcr, respectively. In (b), the ground state is a sin- 
glet on the dotted horizontal line starting from Tr = A^, 
and the transition to the doublet state occurs discontin- 
uously at finite U at the end of the dotted line. These 
examples show that the quantum phase transition can 
be driven by (j) if the other parameters are in the region 
surrounded by the lines for </> = and that for ip — tt. 

We next study how the Coulomb interaction affects the 
Josephson current. To this end, we start with noninter- 
acting case. In Fig. 7 the Josephson current for U = 0, 
which is obtained from Eq. (13), is plotted as a func- 
tion of (j) for several values of Tr. Here the parameters 
are taken to be Ad — |Aj^|, and thus at Tr — |A_r| 
the two couplings of the original model become equal 
Tr = Fi since A^ = F^. When the two couphngs are 
close Tr ~ Fi, the amplitude of the current becomes 
large, and the shape of the dependence deviates from 
a simple sinusoidal form. In the figure the dotted line 
is the result for Tr = F^, and it shows a discontinuous 
jump a,t (f) — ±7r. This singularity is caused by the bound 
state at Fermi level. The NRG results of the Josephson 
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Fig. 6. Phase diagrams of the ground state in the electron-hole 
symmetric case ed = —U/2. Here the impurity SC gap is taken 
to be (a) Arf = |Ah|, and (b) A^ = 2|Aij|. At F^ = 0, the phase 
boundary is given by (/ = 2A£j from the atomic-limit solution. 
For (f) = TV, the bound state comes just on the Fermi level when 
[/ = and Vji = Ajj (= T^). In the case of (b), this happens 
in the region of F^ > Fcr = 1.77064|Afl|, and on the dotted 
horizontal line the ground state is a singlet. 



current are shown in Fig. 8 for several values of C/, where 
the parameters are taken to be Vr = 3.77697|Afl,|, and 
Ad — \A.r\. In this figure the results are plotted only for 
< < TT. The feature of the current at negative ip can 
be seen by rotating the figure around the origin, since 
the current is an odd function of cf) with the period 27r. 
For small C/, the ground state is a singlet for the whole 
range of (j) as seen in the results for U < 20|A/f|, and in 
these cases the amplitude of the current decreases with 
increasing U . We see for U — 25| A/j| that the phase tran- 
sition occurs at ~ 0.767r. Then in the doublet state the 
current flows in the opposite direction from that in the 
singlet state. This change is caused by the spin-flip tun- 
neling through the unscreened local moment. ^^' ^^ There- 
fore, the Josephson current can be used to detect the 
quantum phase transition. 

In quantum dots the coupling to the leads is a tunable 
parameter. We have plotted the current as a function of 
(f) in Fig. 9 for several values of F^ taking U to be (a) 
5.0|Ai{| and (b) 15.0|Afl|. Here the gap is chosen to be 
Ad = \A.r\. Both in (a) and (b), the phase transition 
from the singlet state to doublet state occurs when F^ 
decreases. It can be explained from the fact that the cor- 
relation effects are enhanced for small Tr. However, the 
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Fig. 7. The Josephson current for U = and E^ = as a function 
of (j>. The currents are calculated for several values of F^ from 
Eq. (13). Here A^j {= Fl) is taken to be A^ = |Ai{|, so that F^ 
becomes equal to Fj^ for F^/IA^I = 1.0, where the current for 
(f) = ±7r jumps discontinuously. 




Fig. 8. The Josephson current as a function of </) for various values 
of U. Here the parameters are taken to be F^ = 3.77697Afl, 
Ajj = Ajj, and e^ = —U/2. The ground state is a singlet for U < 
20Afl, while it is a doublet for for U > 30| Ajj|. For U = 25|Ab|, 
the quantum phase transition to the doublet ground state occurs 
at ~ 0.767r, and the direction of the current changes due to 
the spin-flip tunneling. '^^ Note that the current for negative (p is 
symmetric around the origin as that in Fig. 7, since the current 
is an odd function of (j> with the period 2tt. 



amplitude of the current depends on another factor. In 
Fig. 9 (a), the amphtude becomes large when the value 
of Fa approaches Fl (= A^). In this case, the match- 
ing condition at the interface is the dominant factor that 
determines the amplitude of the current. On the other 
hand, the amplitude decreases with Tn in Fig. 9 (b). In 
this case, the ground state is a doublet when F^ becomes 
equal to F^, and thus in the singlet state at F^^ > F^ 
the strong electron correlation dominates the matching 
condition to make the amplitude of the current small. 

4-3 Results away from half-filling 

So far, we have discussed the results at half-filling 
ed = —U/2, where the impurity site is singly occupied 
in average. In that case, the magnetic correlations are 
enhanced, while the charge excitations are suppressed. 
Therefore, when the system goes away from half- filling. 




0.2 0.4 0.6 0.^ 
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Fig. 9. The Josephson current as a function of (^ for various values 
of Tn: the onsite Coulomb repulsion is taken to be (a) U = 
5.0|Aj;| and (b) U = 15.0|Afl|. Other parameters are A^ = 
IAhI, and Ed = -U/2. 



the magnetic doublet ground state should be destabi- 
lized. In quantum dots the impurity potential e^ is the 
parameter that can be controlled by the gate voltages. 
We show in the following the NRG results for the ground- 
state properties in the electron-hole asymmetric case. 

In Fig. 10, one example of the phase diagram is plot- 
ted as a function of U and Tji, where e^ = — 2.5| A^j and 
Ad — |Afl|. The model has the electron- hole symmetry 
a.t U ~ 5.0|Aj^'| and the results coincide with that at the 
same U in Fig. 6(a). As expected, the region where the 
doublet state stabilized becomes narrow away from half- 
filling. Particularly the re-entrant region of the doublet 
state, which is seen at half-filling, almost vanishes in the 
case of Fig. 10. At F/j = 0, all the lines for different 4> 
start from Ua = (e^ + A^)/(— e^), where the level cross- 
ing occurs in the atomic limit. Then, for large U, the 
boundaries become flat. We next consider the e^ depen- 
dence of the ground state. In the following two figures, 
the phase diagrams are plotted in a F^j vs e^ plane for 
two different values oi U: it is chosen to he U = 5.0|Afl;| 
and U = 1.0|Afl| in Figs. 11 and 12, respectively. Here 
the impurity SC gaps is taken to be A^ = |A/j;|. The 
figure 11 is an example of the phase diagram for large 
U, where the impurity is half-filled at e^ = — 2.5|A/i;|, 
and the doublet phase vanishes when ed or ed + U ap- 
proaches the Fermi level at /i = 0. The doublet region 
becomes wider with increasing 0, and this tendency is 
maximized at = tt. For small U , the re-entrant doublet 
region appears as an island in the phase diagram Fig. 12, 
where the line for (j) — tt is shown. The magnetic doublet 
state is most stable at half-filling ed — — 0.5|Afl|. In this 
re-entrant region, the bound state lies close to the Fermi 
level. 
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We have also calculated the Josephson current and the 
average number of electrons in the impurity site (n^) . In 
Fig. 13, the results are plotted as a function of e^ for (a) 
r_R = 1.7|A_r| and (b) r_R — 0.3|A_r|. The parameters are 
taken to be = O.Stt, Ad = \Ar\, and U = 5.0\Ar\: the 
ground state for these parameters are shown in Fig. 11. 
In the case of (a), the ground state is a singlet indepen- 
dent of Ed- The current shows a maximum at half- filling 
ed — — 2.5|Afl|, and (rid) decreases monotonically with 
increasing e^. On the other hand, in the case of (b), the 
quantum phase transition occurs near half-filling, where 
the doublet state is stabilized. The expectation values 
show the singular behavior at the critical point, and 
particularly the Josephson current becomes small and 
changes the direction. 

5. Summary 

In summary, we have studied the ground-state prop- 
erties of a model for the quantum dot embedded in the 
Josephson junction making use of the NRG method. To 
carry out precise calculations in a wide parameter range, 
we have introduced a simplified single-channel model 
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Fig. 13. The Josephson current and the number of the electrons 
in the impurity (n^) as a function of ej for (a) F^ = 1.7|Afl| 
and (b) F^ = O.SlAflJ. The parameters are taken to be </) = O.Stt, 
Arf = |Aii|, and U = 5.0|Ajj|. The ground state is a singlet in 
the case (a), while in (b) the quantum transition occurs and the 
doublet state is stabilized around half-filling. 



that has an extra superconducting gap A^ at the impu- 
rity site. This model captures the essence of the Joseph- 
son and Kondo physics, and can be derived from the 
two-channel model when one of the bulk SC gaps is much 
larger than the others \Al\ ^ |A/j|. This kind of asym- 
metry in the Josephson couplings causes the nontrivial 
changes in the superconducting proximity effects on the 
impurity site, and it affects the low-energy spin dynamics 
of the quantum impurity. 

The numerical results for the simplified model show 
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that the Josephson phase tends to disturb the eonduc- 
tion electron screening of the local moment, and desta- 
bilizes the paramagnetic singlet ground state to drive 
a quantum phase transition to the magnetic doublet 
ground state. At the transition point, the direction of 
the Josephson current changes discontinuously due to 
the spin-flip tunneling. Specifically, when the Josephson 
phase is ~ TT and the other parameters are set so that 
the conductance in a normal limit A^, A^ — > is close to 
the unitary limit value, the bound state comes near the 
Fermi level. This zero mode causes a re-entrant magnetic 
phase transition for the asymmetric Josephson couplings. 
It is known for the Anderson impurity in normal met- 
als that the magnetic transition described by the mean- 
field theory is an artifact due to the approximation, i.e., 
quantum fluctuations stabilize the singlet ground state. 
In the superconducting case, however, the gap changes 
the structures of the low-energy excitations, and makes 
the occurrence of the quantum phase transition possi- 
ble, the phase diagrams for which are provided in the 
present paper. The asymmetries in the Josephson cou- 
plings and the phase difference of the gaps add interest- 
ing varieties in the ground-state properties of the Kondo 
system through the changes in the bound state. 
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Appendix A: Functional method 

The impurity part of the partition function can be ex- 
pressed using the path integral formulation as 



^cff - fvrj^Vrj e- 



S-eff = / / drdr' 
Jo Jo 



(A-1) 






KUr,r') V^ir')vlir') 



U 



'-J\t r£ill{r)^„{T) -ij , (A.2) 

KUry) = ^5]{G°,(ic.„)}-^--"(^--'), (A.3) 

where '^^{t) is a Grassmann number. The information 
about the conduction electrons enters into S'eff through 
the noninteracting Green's function G^i^iiUn)- Thus, if 
the different Hamiltonians give the same G°^(ia;„), they 
also describe the same local properties around the impu- 
rity. ^'^ 



Appendix B: Bogoliubov-de Gennes representa- 
tion 

We provide here the expression of the Green's func- 
tion in terms of the Bogoliubov quasiparticles, which 
will be used in Appendix D to carry out the per- 
turbative calculations. The noninteracting Hamiltonian, 
H° = H° -I- Y1,\=L i? ( ^-^ + ^Ja ) ' <^^^ be diagonahzed 
using the Bogoliubov operator ^na] 



Ho 



Eh Y^ llalOa 



XI 51 Enlialna 



const. 
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(B-1) 

Here the label n = is assigned to the bound state in 
the gap, and that for n > corresponds to the excita- 
tions in the continuum En > min(|Ai|, |A/j|). Then the 
impurity Green's function can be written as. 



GlS^i) = U 
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(B-2) 



Here, Udn is the matrix defined by the Unitary transfor- 
mation 
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(B-4) 



and it has a property U dnU\^ = {\udn\ + \vdn\ )l- 
With these quasiparticle states, the SC correlation for 
noninteracting electrons x^ — (did^)'^ can be expressed 
as 



X'd 



Udov*do { 1 - '2f{Eb) } 

J2 UdnvLi-^ - 2f{En)} 



(B-5) 



n>0 



Note that the first term in the right-hand side, which 
corresponds to the contributions of the bound state, van- 
ishes if the bound state appears at the Fermi energy 
Eb = 0. 

Appendix C: Correlation functions for the zero 
mode Eh — 

As mentioned in §3, the bound state appears just on 
the Fermi level when = tt, F^^ = F^ (= F), and Ed = 0. 



Around the pole at Ei, 
written in the form 



0, the Green's function can be 

(C-l) 
(C-2) 



iO^ 
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r 

|Ai,| 



r 



The spectral weight at the impurity site, a, decreases 
with increasing F because the bound-state wavefunction 
penetrates deep inside the leads for large F. Alterna- 
tively, the Green's function can be expressed using Eq. 
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(B-2) as 



G° 



dd 



(e) 



\Udo\ 



\Vdo\ 
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(C-3) 



Therefore comparing Eqs. (C-2) and (C-3), we have a = 
l^doP + l^'doP- Furthermore, it is concluded that \udo\'^ = 
bdoP = a/2 because of the electron-hole symmetry. Also, 
Udo and Vdo are essentially real at <j) = it. 

In the presence of the zero mode the expression of the 
SC correlation, Eq. (12), can be rewritten in the form 

X'd = 



E- 



Tu:l 



lA, 



/^trdet{G°,(ic.„)}" 
|Al| 



{^l 



lA, 



lAf 



V^?. + |AlP ^^l + \^R? 



(^UIAkP) 



(C-4) 



This expression clearly shows that Xd vanishes for 
I Ail = |Afl|. Note that the behavior of the determi- 
nant near the bound state is deduced form Eq. (C-2) as 

det{G2d(iw„)}"^ ~ -ujl/a^. In Eq. (C-4) this ojI de- 
pendence is canceled out by that in the numerator. 



states is given by 



Es=a — -Es^i 



^^[a^ -2a|x2|] 



=1/2 "" "^ L" ^ ^"|Xd| J • (D-8) 

Thus the level crossing occurs at a — '2\Xd\ ~ ^' where Xd 
is given by Eq. (C-4), and this equation determines the 
phase boundary between the doublet and singlet states. 
The results are shown in Fig. 1 as a function of the ratio 
X = |Afl|/|AL|. At a; = 0, i.e., in the limit of |Al| -^ oo, 
Eq. (C-4) simplifies for T = as 



X"d = 



do; 
— F 

27r 



v/c^2 + |A^|2 + 2rj (v/^2 + |A^|2 + |A^| 

1 
2F^ , (D-9) 



and the critical value of F is given by Fcr/IAi^^j ~ 
1.77064 - - - . In the opposite limit at x — 1 — 0+, i.e., for 
I Ail = |Afl|, the critical value diverges exponentially as 
Fcr/IAfll ~ 4e[V(i-^) + i] ^ so that the doublet ground 
state is stabilized for any values of F. 



Appendix D: Perturbation theory for the zero 
mode 

In the presence of the zero mode, i.e., for (p — tTjTl ~ 
Ffl (= F), and Ed = 0, the following four states are 
degenerate; 



in) 
|iii) 

|IV) 



= |o), 

= 7011 



7^ JO), 
70170x1 



(D-1) 
(D-2) 
(D-3) 
(D-4) 



Here |0) is the vacuum state with respect to the Bo- 
goliubov particles 7„ct|0) = 0. An infinitesimal repulsion 
H^ ~ {U/2) (ud — 1) lifts the degeneracy of the ground 
state. We have calculated the matrix elements {a\H^\a') 
in this subspace making use of Eq. (B-3) and the equa- 
tions provided in Appendix C. We found that the off- 
diagonal elements are zero 



(«|Wd 



= , for a j^ a' 



(D-5) 



and the energy shift is determined simply by the diagonal 
elements 

(I|W^|I) - (II|H^|II) - ^ K - 2ax"dSgniudoVdo)] , 

(D-6) 



(IVIH^IIV) - (II|H^|II) 



U r 2 



2aXdSgn(udoWdo)] • 
(D-7) 

Note that Es=i/2 = (n|7i:^|II) = (III|7i:^|III) is the en- 
ergy shift for the doublet states. Consequently, the en- 
ergy difference between the doublet and the lowest singlet 
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